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ABSTRACT 


An  analytical  method  Is  presented  for  predicting  the  added 
resistance  (thrust  deduction)  arising  from  propeller-hull  in- 
teraction. The  theory  Is  formulated  In  terms  of  the  potential 
flow  about  the  hull  and  appendages  which  are  represented  by 
surface  singularity  distributions.  The  influence  of  the  pro- 
peller is  derived  from  lifting-surface  theory  Including  the 
effects  of  blade  number,  thickness,  skew,  rake,  and  radial  and 
chordwlse  load  distribution.  In  order  to  determine  the  inter- 
action force,  propeller-induced  velocities  and  the  modified  hull 
pressure  distribution  are  computed  with  appropriate  corrections 
to  the  hull  singularities.  The  axial  force  is  then  derived  by 
Integration  of  the  pressure  on  the  hull  and  also  by  application 
of  the  steady-flow  Lagally  theorem. 

The  usefulness  of  this  technique  Is  illustrated  by  its  appli- 
cation to  several  stern  propeller-body-of-revolution  conflgura- 
tlonc.  It  is  shown  that  stern  appendages  contribute  up  to  25 
percent  of  the  total  thrust  deduction.  The  relative  contribu- 
tions of  propeller  loading  and  thickness  are  examined.  It  is 
found  that  the  lifting-surface  representation  predicts  10  to  20 
percent  lower  thrust  deduction  than  the  classical  lifting-line 
(sink  disk)  approximation.  Calculations  for  a series  of  four 
raked  propellers  illustrate  the  significant  (over  50  percent) 
attenuation  of  the  interaction  force  as  rake  is  increased.  It 
is  concluded  that  the  method  is  useful  for  both  the  analysis  of 
a given  design  and  for  parametric  investigations  of  higher 
efficiency  propeller-afterbody  configurations.  The  method  may 
also  be  extended  to  treat  contrarotating  and  ducted  propellers. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  in-house  independent  research  and 
exploratory  development  program  of  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC)  (Work  Unit  No.  1552-119)  and  the  High 
Speed  Submarine  Direct  Laboratory  Funding  Program  (Work  Unit  No.  1520-004 
and  1500-200). 


INTRODUCTION 

The  interaction  force  arising  from  propellers  operating  in  close  proxi- 
mity to  the  ship's  stern  is  a familiar  concept  to  naval  architects.  The 
propeller  accelerates  the  flow  over  the  hull  afterbody.  For  sufficiently 
fine  ship  forms  where  flow  separation  effects  are  minimal,  the  velocity 
increase  and  accompanying  reduction  in  pressure  increases  the  hull  pressure 
drag.  The  higher  velocity  also  increases  the  wall  shear  stress,  and  hence. 
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the  frictional  resistance.  The  net  result  is  that  the  delivered  propeller 
thrust  must  be  greater  than  the  hull  resistance  in  the  absence  of  the 
propeller. 

This  increase  in  resistance  due  to  the  propeller-hull  interaction  is 
defined  in  terms  of  the  thrust  deduction  fraction  t. 


where  is  the  bare  hull  resistance  and  T is  the  propeller  thrust.  The 
thrust  deduction  must  be  known  in  advance  so  that  a propeller  design  will 
meet  the  specified  propulsion  requirements.  One  approach  is  to  conduct 
model-scale  propulsion  tests  using  a stock  propeller  with  similar  principal 
characteristics.  While  this  technique  has  proven  reasonably  satisfactory 
for  many  conventional  designs,  a large  number  of  experiments  are  required 
to  investigate  the  effects  of  different  afterbody  forms,  propeller  loca- 
tions, blade  geometries,  and  loading  characteristics.  Thus,  an  analytical 
prediction  technique  is  desirable  both  from  the  standpoint  of  predicting  the 
interaction  force  for  a given  propeller  and  hull  design,  and  for  economi- 
cally investigating  more  efficient  propeller-hull  configurations. 

Various  techniques  for  the  analysis  and  prediction  of  thrust  deduction 

have  been  reported  in  the  last  forty  years,  as  shown  in  the  comprehensive 
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bibliography  presented  by  Nowackl  and  Sharma.  Dlckman  was  the  first 
investigator  to  provide  a reasonable  theoretical  analysis  of  the  inter- 
action force  between  a hull  and  propeller.  To  represent  an  axisymmetric 
body,  he  applied  the  method  of  discrete  singularities  on  the  body  axis 
together  with  a single  point  sink  to  represent  the  propeller.  With  this 
model,  it  was  possible  to  relate  the  thrust  deduction  to  the  thrust  loading 
coefficient.  During  the  1940's,  some  of  Dickman's  ideas  were  extended,  as 


^Nowackl,  H.  and  S.D.  Sharma,  "Free  Surface  Effects  in  Hull  Propeller 
Interaction,"  The  University  of  Michigan  College  of  Engineering  Report  112 
(Sep  1971).  A complete  listing  of  references  is  given  on  pages  66-67. 

2 

Dickmann,  H.E.,  "The  Interaction  between  Propeller  and  Ship  with  Special 
Consideration  to  the  Influence  of  Waves,"  Jahrbuch  der  Schif fbautechnischen 
Gesselschaft,  Vol.  40  (1939). 
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outlined  by  Weinblura  in  his  survey  paper.  Korvin-Kroukovsky  also  used  the 
method  of  singularities,  but  used  a constant-strength  sink  disk  as  a pro- 
peller model,  and  corrected  the  body  sources  to  account  for  both  the  pro- 
peller induced  flow  and  the  boundary  layer  displacement  thickness. 

More  recent  developments  have  logically  followed  from  advances  In  pro- 
peller theory  and  the  advent  of  high-speed  digital  computing  capabilities. 
Thus,  Beveridge^  was  the  first  to  apply  the  Douglas-Neumann  three-dimension- 
al potential  flow  calculation  to  represent  the  body.  In  a later  investiga- 
tion, Beveridge^  introduced  a sink  disk  with  radially  varying  strength 
derived  from  propeller  lifting-line  theory.  Based  on  this  method,  predic- 
tions for  three  different  realistic  hull  forms.  Including  one  surface  ship, 

compared  reasonably  well  with  experimental  data.  A similar  technique  has 

7 8 

also  been  applied  to  contrarotating  propellers.  ’ 

In  the  present  work,  an  extended  treatment  of  the  potential  flow  analy- 
sis is  developed  which  differs  in  two  fundamental  respects  from  previous  ap- 
proaches. First,  a more  comprehensive  and  realistic  representation  of  the 
propeller  is  introduced,  based  on  a lifting-surface  formulation.  This  over- 
comes the  limitations  of  the  lifting-line  sink-disk  approximation  by 
Including  the  additional  effects  of  blade  thickness,  skew,  rake,  and  chord- 
wise  load  distribution.  As  such,  the  propeller  calculation  is  comparable  in 
scope  and  accuracy  to  currently  available  lifting-surface  design  methods. 
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Weinblum,  G.,  "The  Thrust  Deduction,"  American  Society  of  Naval  Engin- 
eers, Vol.  63  (1951). 

^Korvin-Kroukovsky,  B.V. , "Stern  Propeller  Interaction  with  a Stream- 
line Body  of  Revolution,"  International  Shipbuilding  Progress,  Vol.  3,  No. 
17  (1956). 

^Beveridge,  J.L.,  "Pressure  Distribution  on  Towed  and  Propelled  Stream- 
line Bodies  of  Revolution  at  Deep  Submergence,"  David  Taylor  Model  Basin 
Report  1665  (Jun  1966). 

^Beveridge,  J.L.,  "Analytical  Prediction  of  Thrust  Deduction  for  Sub- 
mersibles  and  Surface  Ships,"  Journal  of  Ship  Research,  Vol.  13,  No.  4 (Dec 
1969). 

^Nelson,  D.M. , "Development  and  Application  of  a Lifting-Surface  Design 
Method  for  Counterrotating  Propellers,"  Naval  Undersea  Center  TP  326  (Nov 
1972). 

g 

Beveridge,  J.L.,  "Thrust  Deduction  in  Contrarotating  Propellers,"  Naval 
Ship  Research  and  Development  Center  Report  4332  (Nov  1974). 
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In  previous  Investigations,  the  thrust  deduction  has  been  derived  by 
applying  the  Legally  steady-flow  theorem  to  the  propeller  sink-disk  singu- 
larities. In  the  present  formulation,  It  is  more  convenient  to  consider 
the  body  flow  directly.  Thus,  time-averaged  propeller  induced  velocities 
and  modified  hull  pressure  distributions  are  computed,  including  appropriate 
corrections  to  the  body  singularity  strengths.  The  force  Is  then  calculated 
by  integrating  the  pressure  and  also  by  applying  Lagally's  theorem  to  the 
body  singularities.  It  Is  assumed,  as  before,  that  for  a given  representa- 
tion of  the  propeller,  the  added  hull  resistance  arises  entirely  from  the 
reduction  In  afterbody  pressure  and  that  this  pressure  distribution  may  be 
derived  solely  from  potential  flow  considerations.  Although  it  is  recog- 
nized that  the  boundary  layer  at  the  stern  is  relatively  thick,  the  poten- 
tial flow  formulation  has  been  widely  accepted  on  the  basis  of  agreement 
between  predicted  and  measured  thrust  deduction.  Recently,  wind-tunnel 

experiments  were  conducted  on  streamlined  bodies  of  revolution  with  and 

9 

without  a propeller  in  operation.  The  results  show  that  while  the  theory 
cannot  satisfactorily  predict  the  absolute  pressure  distributions  near  the 
stern,  the  difference  in  pressure  due  to  the  action  of  the  propeller  is 
predicted  remarkably  well.  It  was  also  found  that  increases  in  wall  shear 
stress  contribute  less  than  5 percent  of  the  integrated  pressure  force. 

For  these  reasons,  it  appears  that  the  Douglas-Neumann  potential-flow  calcu- 
lation is  a sound  approach,  at  least  for  nonseparating  hull  forms.  Moreover, 
the  calculation  of  the  detailed  pressure  distribution  will  serve  as  a 
necessary  first  step  in  future  treatments  of  the  viscous  flow  problem. 

In  this  report,  the  theoretical  basis  and  numerical  techniques  for 
predicting  the  thrust  deduction  are  presented.  The  analysis  is  restricted 
to  deeply  submerged  bodies,  for  which  the  hull  potential  flow  calculation 
is  only  briefly  reviewed,  being  extensively  documented  in  the  cited 
literature.  (The  theory  can,  in  principle,  be  extended  to  surface  ship 
applications;  the  free-surface  would  be  approximately  represented  by 


q 

Huang,  T.  et  al.,  "Propeller/Stern/Boundary-Layer  Interaction  on 
Axisymmetric  Bodies:  Theory  and  Experiment,"  DTNSRDC  Report  76-0113 
(Dec  1976). 


4 


reflecting  the  hull  and  propeller  (images)  in  the  waterline  plane.)  The 
propeller  representation  and  field  point  velocity  calculations  are  dis- 
cussed in  some  detail.  The  further  analysis  considers  the  determination  of 
the  modified  hull  pressure  and  solution  for  the  interaction  force  by 
pressure  Integration  and  by  application  of  Lagally's  theorem. 

Comparisons  between  theoretical  and  experimental  thrust  deduction  are 
given  for  deeply  submerged  stern  propeller  body-of-revolution  configura- 
tions both  with  and  without  stern  appendages.  These  configurations  were 
chosen  for  initial  calculations  because  the  body  geometry  characteristics 
afford  easier  computation.  However,  it  is  noted  that  in  application  to 
submarines  and  torpedoes,  the  thrust  deduction  is  of  great  practical  im- 
portance in  selecting  propeller  characteristics  (e.g.,  diameter)  for  maxi- 
mum propulsive  efficiency.  In  the  examples  presented,  it  is  shown  that 
stern  appendages  develop  as  much  as  25  percent  of  the  thrust  deduction. 

The  relative  contributions  of  propeller  loading  and  thickness  are  examined 
and  compared  with  the  classical  lifting-line  sink-disk  results.  It  is 
found  that  lifting-surface  effects  reduce  the  thrust  deduction  by  10  to  20 
percent.  Calculations  for  a series  of  four  rakes  illustrate  the  significant 
(over  50  percent)  attenuation  of  the  interaction  force  as  rake  is  Increased. 
Based  on  these  examples,  it  is  concluded  that  the  theory  provides  a useful 
technique  for  both  the  analysis  of  a given  design  and  for  parametric 
studies  of  higher  efficiency  propeller-hull  configurations. 

PROPELLER-HULL  INTERACTION  IN  POTENTIAL  FLOW 

SINGULARITIES  REPRESENTING  THE  HULL  AND  APPENDAGES 

The  propeller-hull  interaction  analysis  rests  on  the  computation  of 
the  potential  flow  about  the  hull  in  the  presence  of  the  propeller.  A 
great  deal  of  effort  has  been  devoted  in  the  past  to  developing  accurate 
and  computationally  efficient  techniques  for  calculating  the  flow  about 


5 


r 


arbitrary  three-dimensional  bodies.  In  the  method  of  Hess  and  Smith,  ’ 

the  body  surface  is  approximated  by  planar  quadrilateral  elements  and  the 

solution  is  derived  in  terms  of  simple  source  distributions.  A recently 
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developed  computer  code  based  on  this  approach  is  employed  in  the  present 
work.  Briefly,  the  formulation  is  as  follows. 

The  body  Is  assumed  to  be  deeply  submerged  and  advancing  at  a constant 
velocity,  V^,  in  an  incompressible,  Invlscid  fluid.  By  considering  only  the 
time-averaged  propeller  disturbance  field,  the  flow  is  steady  relative  to  a 
Cartesian  coordinate  system  r = (x,y,z)  advancing  with  the  body.  Outside  the 
propeller  blade  row  and  slipstream  the  flow  is  irrotational.  Thus,  a velo- 
city potential  (Ji(r)  exists  such  that  V(r^)  » V((>(r)  and 

^(r)  = 0 (1) 

where  “r  is  outside  the  body  and  the  propeller  slipstream.  The  boundary 
conditions  are  that  the  velocity  must  be  tangent  to  the  body  surface, 

n(rg)  • V(t>(rg)  = 0 (2) 

and  approach  the  free  stream  value  at  a large  distance  from  the  body- 
propeller  system 

V(|)(^  V^,  as  l‘r|  -»■  “ (3) 


A solution  for  (p(r)  which  satisfies  equations  (1)  and  (3)  may  be  written  in 
terms  of  a surface  source  distribution,  a(r„),  as 

D 


Mr)  = - 


4it 


r - r. 


dS(rg)  + r • 


(4) 


^^Hess,  J.L.  and  A.M.O.  Smith,  "Calculation  of  Nonlifting  Potential  Flow 
About  Arbitrary  Three-Dimensional  Bodies,"  Journal  of  Ship  Research  (Sep 
1964). 

^^Hess,  J.L.  and  A.M.O.  Smith,  "Calculation  of  Potential  Flow  About  Arbi- 
trary Bodies,"  Pergamon  Press,  Progress  in  Aeronautical  Sciences,  Vol.  8 
(1966). 

^^Dawson,  C.W.  and  J.S.  Dean,  "The  XYZ  Potential  Flow  Program,"  NSRDC 
Report  3892  (Jun  1972). 
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Here  't’pCr)  is  the  potential  due  to  the  propeller  which  satisfies  the 
condition 


V(J.p(r)  0 lr|  -»■  “ (5) 

Inserting  this  expression  for  iJ'Cr)  into  equation  (2)  yields  an  Integral 
equation  for  the  unknown  source  strengths 

h ^ IV’^bI 


It  may  be  observed  that  the  propeller  presents  a modified  onset  flow  to  the 
body  (right-hand  side  of  equation  (6)).  Moreover,  the  change  in  source 
strength  caused  by  the  presence  of  the  propeller  depends  only  on  the  com- 
ponent of  Induced  velocity  normal  to  the  body  surface. 

A numerical  solution  of  equation  (6)  is  obtained  by  representing  the 
body  surface  using  planar  quadrilateral  elements.  It  is  assumed  that  the 
source  density  is  constant  over  each  element  and  the  integral  equation  is 
replaced  by  a set  of  linear  algebraic  equations. 


E 


0 . 
1 


ij 


where  the  coefficients. 


are  given  by 


C.  . 
ij 


<5,.  r n.’V  ,3 — — I dS. 

i J 


element 


(7) 


(8) 


and  is  the  onset  flow  evaluated  at  a selected  control  point  (e.g.,  cen- 
troid) of  each  quadrilateral 


V 


1 


[V„  + 


(9) 


Since  the  coefficients  C..  depend  only  on  the  body  geometry, 
-1 

matrix  C , need  only  be  derived  once  for  a given  hull  form. 


the  inverse 
It  is  then 
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possible  to  rapidly  compute  the  source  strengths  corresponding  to  a 
number  of  propeller  onset  flows  as 


<^1  = E c 


V 

Ij  i 


(10) 


Once  the  body  source  strengths  are  known,  it  is  straightforward  to 
compute  velocities  and  pressures  at  points  on  the  body  surface  or  in  the 
surrounding  field.  It  is  also  possible  to  determine  the  resultant  force 
acting  on  the  body  as  described  below. 

SOLUTION  FOR  THE  THRUST  DEDUCTION 

In  general,  the  influence  of  a stern  propeller  causes  a net  increase 
in  the  hull  resistance.  Two  methods  are  available  to  compute  the  force 
exerted  on  the  body.  In  the  first  approach,  the  axial  force,  F^,  is 
derived  by  integrating  the  pressure  over  the  body  surface. 


F = e 

X X 


•i 


pCV  - 

with 

propeller 


p(^b) 

without 

propeller 


n(fg)  dS(fg) 


(11) 


where  the  pressure  is  found  from  Bernoulli's  equation 


p = --^V  • V 


(12) 


The  velocity  is  computed  at  the  control  points  of  each  quadrilateral  as 


I"b.-"b, 
1 j 


AS.  + V + V(})  (r„  ) 
j P B,. 


A distinct  advantage  of  this  method  is  that  the  effect  of  body  form  can  be 
determined  by  examining  the  distribution  of  the  pressure  integral.  Also,  a 
detailed  knowledge  of  the  pressure  distribution  is  an  Important  first  step 
in  solving  the  viscous  flow  over  the  afterbody. 
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The  interaction  force  may  also  be  derived  by  application  of  Lagally's 
13  14 

theory  ’ for  a body  immersed  in  a steady  potential  flow.  The  solution 
for  the  force  on  a body  due  to  an  isolated  point  source  is  developed  in 
Appendix  A as 


■ 


m D 


(13) 


where  is  the  undisturbed  onset  velocity  generated  by  the  singularity. 

This  equation  is  also  valid  for  a point  doublet  singularity.  It  will  be 
shown  subsequently  that  the  propeller  disturbance  arises  from  suitable 
distributions  of  sources  (blade  thickness)  and  doublets  or  equivalently, 
line  vortices  (blade  loading).  Thus,  the  axial  force  arising  from  the 
propeller-hull  interaction  may  be  written  as 


F 

X 


e 

X 


) V<^  (^  ) dS( 
P “ 


or  in  discrete  form 


(14) 


F 

X 


E 


(15) 


If  only  the  total  force  is  required,  this  form  is  simpler  for  computation 
and,  in  any  case,  provides  a convenient  check  in  the  numerical  evaluation 
of  equation  (11). 

Once  the  interaction  force  is  found,  the  thrust  deduction  fraction,  t, 
is  given  by 

F 

t = ^ (16) 


13 

Cummins,  W.E.,  "The  Force  and  Moment  on  a Body  in 
Potential  Flow,"  Journal  of  Ship  Research,  Vol.  1,  No. 

^^Milne-Thomson,  L.M. , "Theoretical  Hydrodynamics," 
New  York,  N.Y.,  2nd  edition  (1950). 
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The  Macmillan  Company, 
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where  T is  the  propeller  thrust.  Up  to  this  point,  it  has  been  assumed 
that  the  propeller  has  been  represented  by  an  appropriate  distribution  of 
singularities  external  to  the  body.  Within  the  framework  of  the  potential 
flow  formulation,  the  solution  to  the  interaction  problem  is  derived  com- 
pletely in  terms  of  propeller  disturbance  velocities  at  points  on  the  body 
surface.  It  is  now  appropriate  to  set  forth  the  theoretical  basis  and 
numerical  techniques  for  calculating  these  velocities. 

PROPELLER  FIELD  POINT  VELOCITIES 

In  the  foregoing  analysis,  the  modified  body  flow  in  the  presence  of  a 
propeller  is  derived  in  terms  of  induced  velocities  on  the  body  boundary. 

It  is  primarily  in  the  treatment  of  the  propeller  that  the  present  inter- 
action analysis  differs  from  earlier  investigations.  Owing  largely  to 
advances  in  design  theory  and  high-speed  computing  capabilities,  it  has 
been  possible  to  introduce  a more  realistic  analytical  representation  of 
the  propeller. 

Previously,  the  propeller  was  approximated  as  a sink  disk.  In  that 
model,  the  diameter,  axial  location,  and  radial  distribution  of  loading  are 
explicitly  represented.  In  fact,  as  will  be  shown  shortly,  a sink  disk  of 
strength 


, . _ Z dr  1 , . 

ZirJ  dr'  r'  tan  6.  (r"*)  ^ 

r 1 

generates  the  circumferential  average  velocity  field  of  a moderately  loaded 
lifting-line  representation  of  a propeller  with  bound  circulation  F,  pitch 
Ztt  r tan  3^  (r) , and  Z symmetrically  spaced  blades.  Beveridge  has  demon- 
strated that  this  model  satisfactorily  predicts  the  thrust  deduction  for 
conventional  propeller  geometries.  However,  it  is  evident  that  this  sim- 
plification breaks  down  for  raked  propellers  in  which  the  blade  sections 
are  displaced  axially.  Moreover,  the  effects  of  blade  thickness  and  finite 
chordlength,  while  perhaps  of  small  consequence  to  the  thrust  deduction, 
may  be  Important  to  the  body  pressure  distribution  and  boundary  layer 
characteristics  in  the  immediate  vicinity  of  the  propeller.  In  view  of 
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these  limitations,  a more  complete  representation  of  the  propeller  based 
on  lifting  surface  theory  is  required. 

A number  of  propeller  lifting  surface  theories  have  been  developed 
and  programmed  for  computer-aided  design  calculations.  The  method  of 
Kerwin^^’^^  is  representative  of  the  state-of-the-art  and  was  selected  for 
the  present  application  because  the  theory  and  numerical  techniques  have 
been  extended  to  calculate  Induced  velocities  at  arbitrary  field  points. 
Free  space  pressure  predictions  based  on  this  method  show  excellent  agree- 
ment with  experimental  measurements.^^  Both  the  theory  and  numerical 
analysis  are  conveniently  divided  into  a lifting-line  analysis  and  lifting- 
surface  corrections  arising  from  blade  thickness,  blade  location  (skew  and 
rake),  and  chordwise  variation  in  loading.  Therefore,  it  is  possible  to 
examine  the  separate  contributions  of  various  propeller  characteristics  to 
the  thrust  deduction. 

LIFTING-LINE  THEORY 

The  basis  of  analytical  propeller  design  methods  is  the  moderately 

18  19 

loaded  lifting  line  theory.  ’ The  analysis  considers  the  flow  field 
associated  with  the  steady  loading  on  a propeller  with  symmetrically  spaced 
blades.  In  accordance  with  circulation  theory,  the  pressure  loading  on 
the  blades  arising  from  camber  and  incidence  can  be  represented  by 
distributions  of  bound  and  free  vorticity.  In  the  lifting  line  approxima- 
tion, each  blade  is  replaced  by  a single  concentrated  vortex  line  with 


^^Kerwin,  J.E.  and  R.  Leopold,  "A  Design  Theory  for  Subcavitatlng 
Propellers,"  Transactions  SNAME,  Vol.  72  (1964). 

^^Kerwln,  J.E. , "Computer  Technique  for  Propeller  Blade  Section  Design," 
International  Shipbuilding  Progress,  Vol.  20,  No.  227  (Jul  1973). 

^^Denny,  S.B.,  "Comparisons  of  Experimentally  Determined  and  Theoreti- 
cally Predicted  Pressures  in  the  Vicinity  of  a Marine  Propeller,"  Naval 
Ship  Research  and  Development  Center  Report  2349  (May  1967). 

18 

Lerbs,  H.W. , "Moderately  Loaded  Propellers  with  Finite  Numbers  of 
Blades  and  an  Arbitrary  Distribution  of  Circulation,"  Transactions  SNAME, 
Vol.  60  (1952). 

19 

Morgan,  W.B.  and  J.W.  Wrench,  "Some  Computational  Aspects  of  Propeller 
Design,"  Methods  in  Computational  Physics,  Vol.  4,  Academic  Press  Inc.,  New 
York,  N.Y.  (1965). 
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bound  circulation,  r(r).  Conservation  of  vortlcity  requires  that  a free- 

d r c ) 

vortex  sheet  of  strength  - - ~ shed  from  each  lifting  line.  These 
sheets  are  assumed  to  lie  on  the  surfaces 


x(r,>fi)  = r tan  B^(r)  Rjj  1 r <.  R 

k - 1,2, ...Z 

where  is  a cylindrical  coordinate  system  fixed  to  the  propeller  as 

shown  in  Figure  1.  angular  position  of  the  blade. 


k Z 


1,2, . . .Z 


and  3^(r)  is  the  flow  angle  defined  in  terms  of  the  relative  velocity  at 
the  blade  section  as  shown  in  Figure  2. 


tan  B^(r) 


Va(r)  + u'(0,r,y>j^) 
nr  - u^(0,r,<Fj^) 


(20) 


Expressions  for  the  induced  velocity  u'(x,r,<^)  = (u',u',u^)  can  be 
derived  by  applying  the  Biot-Savart  law  to  the  distributions  of  vortici- 
ty.^^’^^  The  Induced  velocity  due  to  the  bound  vortex  lines,  u^,  is 


where 

(0,  cos  sin 

(x,  r cos  - p cos  r sin  <fi  - p sdn 


(21) 


^®Wu,  T.Y.,  "Some  Recent  Developments  in  Propeller  Theory,"  Schiffs- 
technik,  Vol.  9,  No.  47  (1962). 
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Figure  1 - Coordinate  Systems  for  Propeller  Lifting  Line  Theory 


r 


Similarly,  the  velocity  induced  by  the  free-vortex  sheets,  uf,  is  given  by 


e X D 


da  dp  (22) 
J 


k 's' 


where 


(X^(p),  - p cos  a,  p sin  u) 


UJ(p)  + p^ 


= (x  - X^(p)  a,  r cos  <fi  - p cos  a,  r sin  V’  - p sin  a) 


X^(p)  = p tan  0^(p) 


If  these  expressions  are  evaluated  at  the  lifting  line  (0,  r,  and 
inserted  into  equation  (20),  an  integral  equation  is  obtained  relating 
r(r)  and  6^(r). 

In  design  applications  where  a prescribed  thrust,  T,  is  to  be  developed, 
a second  relationship  can  be  derived  by  applying  the  Kutta-Joukowski  law 
(with  an  empirical  correction  for  profile  drag  - (see  Figure  2)  yielding 


T = Z 


= pZ  J r(r)  [fir  - u^(0,r,V>^)]  [l  - e(r)  tan  B^(r)]  dr 


The  solution  is  normally  found  by  an  iterative  procedure,  starting  with  an 
estimate  of  tan  r(r)  is  then  computed  from  equations  (20)  and  (22)  and 

used  to  calculate  a thrust  according  to  equation  (23).  This  procedure  is 
repeated  until  the  desired  value  of  thrust  is  obtained.  Highly  efficient 
numerical  techniques  based  on  asymptotic  series  expansions  have  been 
developed  to  perform  these  computations.  In  the  present  work,  it  is  assumed 
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that  both  r(r)  and  B.(r)  have  been  determined,  either  from  design  calcula- 

^ 21  22 
tlons  or  from  a separate  performance  prediction.  ’ 

The  Induced  velocity  varies  with  angular  position,  <P , corresponding  to 

a time  dependence  In  a reference  frame  fixed  to  the  hull,  l.e., 

* n 3t 


By  virtue  of  the  propeller  symmetry,  the  Induced  velocity  may  be  resolved 
into  a time-averaged  or  steady  velocity  ■u(x,r)  and  harmonics  in  blade 
passage  frequency,  Zil.  The  steady  component  or  zero  harmonic  which  gives 
rise  to  the  thrust  deduction,  is  given  by 


u(x,r) 


u'(x,r,V’) 


d'fi 


(24) 


The  separate  contributions  from  the  bound  and  trailing  vorticity  can  be 
determined  by  introducing  expressions  of  the  type 


oo 

= S e cos  n('^  “ ) 

|D  J n=0  " 


f J (kr)  J (kp)e 
Jq  n n 


■k  ixl  "n  “ 1 “ 0 

- 2 n > 0 


and  performing  the  integration  in  equation  (24) . The  following  results  are 
obtained : 
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Cummings,  D.E.,  "Numerical  Prediction  of  Propeller  Characteristics," 
Journal  of  Ship  Research,  Vol.  17,  No.  1 (Mar  1973). 
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Tsao,  S.S.,  "Documentation  of 
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Ug  (x,r) 


Ug  (x,r) 
r 


U (x,r) 


(x,r) 


(x,r) 

r 


^ X “ 

■Si  ^ 


dk  dp  X > 0 


1^1  ji  j 

e dk  dp 


■ ■ l^r  i'  f X ■*''  ■'o<'“>  " 5 » 


u (x.r) 

V 


'R„  •'0 

n 


The  total  tangential  velocity  component  is  given  by 


° 


/ N zr(r)  „ „ 

u/x.r)  = - -y--  X > 0.  Rg  < r < R 
= 0 elsewhere 

It  follows  that  the  steady  velocity  induced  by  a moderately  loaded  lifting 
line  propeller  at  points  outside  the  slipstream  is  given  by 


•^R,,  Jq 


dk  JpCkr)  Jj^(kp)  e dk  dp  x > 0 


-k  X 

e ' ' dk  dp 


u^(x,r)  = 0 


f 


and  arises  solely  from  the  trailing-vortex  sheets.  As  shown  In  Appendix  B, 
this  is  also  the  velocity  field  due  to  a sink  disk  of  strength 


o(r) 


^ r‘'  dr  _1__ 

2t\  dp  X^(p) 


dp 


-2u^(0,r) 


(27) 


The  method  of  Kerwin,  originally  developed  to  compute  all  harmonics 
of  the  induced  velocity,  is  based  on  a direct  numerical  integration  of 
equations  (21)  and  (22).  The  continuous-vortex  sheet  is  divided  into  a set 
of  discrete  helical  vortex  lines  of  constant  strength.  For  computational 
efficiency,  the  integration  interval  is  divided  into  segments  of  increasing 
size  with  distance  from  the  field  point  and  evaluated  using  a five-point 
Gauss  integration  formula.  The  velocity  is  evaluated  at  a selected  number 
of  angular  positions  between  two  blades  and  resolved  into  blade  frequency 
harmonics  by  Fourier  analysis.  Only  the  zero  harmonic  is  used  in  the 
present  application.  As  an  example  of  the  numerical  accuracy  of  this  method, 
a comparison  between  the  computed  and  exact  solutions  of  the  axial  velocity 
at  the  propeller  disk  plane  for  two  selected  loading  distributions  are  shown 
in  Figure  3.  The  velocity  computed  at  the  lifting  line  is  also  shown  for 
comparison  with  the  circumferential  average. 

At  distances  greater  than  about  one  radius,  the  propeller  disturbance 
velocity  is  generated  essentially  by  the  lifting-line  sink  disk.  In  order 
to  more  accurately  derive  the  near-field  influence  of  the  propeller,  the 
lifting-surface  representation  must  be  applied. 


LIFTING-SURFACE  CORRECTIONS 

Propeller  lifting-surface  analyses  have  been  developed  as  a logical 
extension  of  thin  planar  hydrofoil  theory.  Thus,  blade  loading  and  thick- 
ness are  represented  in  terms  of  vorticity  and  source-sink  distributions 
located  on  a reference  surface  approximating  the  actual  blade.  This  is 
normally  taken  to  be  the  pitch  surface  derived  from  lifting-line  calcula- 
tions. Within  the  blade  outline  a distribution  of  bound  and  free  vorticity 
is  established  to  represent  both  the  radial  and  chordwise  load  variation. 

A free-vortex  sheet,  previously  assumed  to  be  shed  from  the  lifting  line, 
now  originates  from  the  blade  trailing  edge. 
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Figure  3 - Field  Point  Velocity  Calculations  in  Propeller  Plane 
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Figure  3a  - Comparison  of  Field  Point  Velocity  Computations  With 
Exact  Solution  for  ^ tan  $. (r)  = 0.4234  (Constant  Pitch) 
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It  is  convenient  to  consider  the  induced  velocities  as  the  sum  of  the 

lifting-line  result  and  separate  lifting-surface  corrections.  The  induced 

velocity  due  to  the  bound  and  free  vorticity,  "e  Y„(r,v>)  and  "e  Y (r,^),  may 

r D s s 

be  written  as 

^ (p) 

e ^ D 

YR(PfOi)  — j ^ dp  da  (28) 

“ l«sl 


and 


u,'.(x,r,<,9)  = 


Z 

E 

k=l 


(p)+p 


'fij  (p) 


e X D 


(29) 


where  'P  ,T,  (r)  are  the  blade  leading  and  trailing  edge  positions.  Continu- 
^k  ^ 

ity  of  vorticity  requires  that 


Y (r,<p)  = - ^ f YR(r,a)  J X^(p)+p^  da 

® -Vrr)  ® ^ ^ 


dr 


(30) 


and  hence  u',  can  be  rewritten  as 


u'(x,r,'P)  = E / 


dr(r)  f 

dr  I 


da 


e X D 
s s 


'^T,(P) 

k 


(P) 


Y*(p ,a) 


e X D 
s s 


(31) 
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where 


Y*(r,<^)  - y iT,<p) 
s s 


<f>.  < <(>  < <P. 


The  induced  velocity  due  to  a distribution  of  sources  over  the  blade 


surfaces  is  given  by 


u;(x.r.V)  - - ^ -L  f f 


•^k 


a(p,a)  V dp  da  (33) 

1^  I 


where  a(r,<p)  is  the  source  strength  density. 

It  follows  that  the  lifting  surface  induced  velocity  may  be  considered 


as  the  sum  of  four  terms 


u(x,r,<o)  = u'Cx.r.V)  + 


k=l  •'r,.  ^ 


V^(P)  ^ 

r f k e X D 

(P)+P  J Y*(p,a)  “ 3® 
> (p)  Id  1 


+ [u'(x.r.</>)]^^j^^^g  line  ^ 


where 


© » velocity  due  to  a distribution  of  bound  vorticity  on  the  blade 

surface  (equation  (28)) 

(2)  = velocity  due  to  a distribution  of  free  vorticity  on  the  blade 

surface 

= velocity  due  to  the  lifting  line  free-vortex  sheet  (equation  (12)) 
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= velocity  due  to  a distribution  of  sources 

In  the  method  of  Kerwin^^  (available  as  a FORTRAN  computer  program 
FPV),  the  numerical  evaluation  is  again  based  on  discrete  singularities. 

On  the  blade  surface,  a grid  of  radial  and  helical  lines  is  constructed  to 
form  a lattice  of  source  and  vortex  elements  as  illustrated  in  Figure  4. 

The  elemental  singularity  strengths  are  determined  as  follows: 

1.  Radial  vortex  elements  are  required  to  produce  the  desired  chord- 
wise  load  distribution  and  the  known  bound  circulation  r(r)  at  each  radius. 

2.  Helical  vortex  elements  must  satisfy  conservation  of  vorticity. 

3.  Source  and  sink  elements  are  required  to  generate  the  same  chord- 
wise  velocity  distributions  as  the  section  thickness  form  would  produce  in 
two-dimensional  flow.  Also  the  sum  of  sources  and  sinks  must  equal  zero. 


As  in  the  lifting-line  analysis,  the  velocity  is  computed  at  a set  of  angu- 
lar positions  for  each  field  point  (x,r)  and  resolved  into  blade  frequency 
harmonics. 

It  may  be  noted  that  radial  vortex  lines  (bound  circulation)  do  not 
contribute  to  the  steady  axial  and  radial  velocities  and  need  only  be 
computed  to  determine  the  strength  of  the  helical  (free)  vortices  on  the 
blade.  It  is  also  straightforward  to  show  that  the  steady  velocity  is 


independent  of  the  skew  angle. 


as  would  be  expected  physically. 


On  the  other  hand,  blade  rake,  the  axial  displacement  of  blade  sections  aft 
of  the  lifting  line  plane,  is  of  marked  importance.  This  is  manifested  in 
the  free  vorticity  term  (2)  in  equation  (34)  which  "corrects"  for  the 
starting  position  of  the  traillng-vortex  sheets.  These  and  other  features 
of  the  propeller  calculations  are  best  illustrated  by  considering  specific 
examples. 


EXAMPLE  CALCULATIONS  AND  COMPARISON  WITH  EXPERIMENTS 
The  equations  derived  in  the  foregoing  theoretical  analyses  have  been 
programmed  for  computer-aided  numerical  solution.  The  calculation  is  per- 
formed by  interfacing  three  separate  programs:  (1)  the  hull  potential  flow 
solution  (PF) , (2)  propeller  field  point  velocities  (FPV),  and  (3)  the 
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1 


SOURCE  AND  VORTEX  LINE  MODEL  OF  PROPELLER  BLADE 


Figure  4 - Representation  of  Propeller  Blade 
Line  Lines  (Projected  View  Looking 


By  Source  and  Vortex 
Forward) 
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interaction  analysis  (CALCTD) • A block  diagram  illustrating  this  procedure 
is  given  in  Figure  5.  It  should  be  noted  that  the  first,  and  most  time- 
consuming  task,  is  to  assemble  the  necessary  hull  offsets  and  propeller 
geometry  and  loading  data  in  a suitable  form. 

Example  calculations  have  been  conducted  for  several  propeller  body- 
of-revolutlon  configurations,  with  and  without  cruciform  stern  appendages. 
These  examples  were  chosen  to  illustrate  important  features  of  the  analy- 
tical method  and  to  provide  experimental  verification  of  the  theory.  A 
summary  of  the  results  is  presented  in  Table  1 showing  pertinent  character- 
istics of  each  propeller-hull  configuration  and  a comparison  of  predicted 
and  measured  thrust  deduction  fractions.  Detailed  numerical  results  and 
discussions  of  each  example  are  given  in  the  following  sections. 

EXAMPLE  1:  APPENDED  SERIES  58  FORM 

As  a first  check  on  the  computational  procedure,  a Series  58  Form 
originally  calculated  by  Beveridge^  was  selected  for  analysis.  The  hull  is 
a streamlined  body  of  revolution,  DTNSRDC  Model  4620,  fitted  with  cruciform 
stern  appendages.  The  propeller,  DTNSRDC  3638,  is  a 5-bladed  wake-adapted 
design  located  98  percent  of  the  hull  length  from  the  bow.  Offsets  and 
particulars  of  the  hull  are  listed  in  Table  2.  A drawing  of  the  propeller 
is  given  in  Figure  6. 

The  quadrilateral  representation  of  the  hull  and  appendages  is  illus- 
trated in  Figure  7 (Identical  to  that  used  by  Beveridge) . Note  that  for 
ease  of  computation,  the  horizontal  control  surfaces  are  also  used  to 
represent  the  upper  and  lower  rudders  and  the  forebody  is  approximated  by 
reflecting  the  afterbody  about  the  hull  midlength  (this  latter  simplifica- 
tion is  shown  to  be  valid  in  later  examples).  The  propeller  circulation 
and  hydrodynamic  pitch  distributions,  r(r)  and  2Trr  tan  6j|^(r),  were  obtained 
from  lifting-line  calculations  and  are  shown  in  Figure  8.  The  axial  and 
radial  propeller  field  point  velocities  induced  at  control  points  on  the 
body  boundary  are  presented  in  Figure  9.  Since  the  propeller  has  no  rake, 
the  lifting-surface  corrections  arise  solely  from  blade  thickness  and  chord- 
wise  load  distribution.  These  effects  are  seen  to  decay  more  rapidly  than 
the  lifting-line  disturbance  field,  contributing  less  than  10  percent  at 
distances  beyond  one  propeller  radius. 
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Figure  5 - Interaction  Analysis  - Computational  Procedure 
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TABLE  2 - OFFSETS  AND  PARTICULARS  FOR  SERIES  58  FORM,  MODEL  4620^’^^ 


X 

Y 

X 

— 

Y 

X 

in. 

y 

in. 

X 

in. 

y 

in. 

0,00 

000.0 

0.0000 

0.000 

0.52 

0.4818 

11.82 

3.6 

0, 1427 

3.500 

0.54 

0.4755 

11.66 

7.2 

0.2029 

4.977 

0.56 

100.8 

0.4684 

11.49 

0.06 

10.8 

0.2490 

6.108 

0.58 

104.4 

0.4603 

11.29 

0.08 

14.4 

0.2873 

7.047 

0 60 

108.0 

0.4513 

11.07 

0. 10 

18.0 

0.3200 

7.850 

0.62 

111.6 

0.4414 

10.83 

0. 12 

21,6 

0.3485 

8.549 

0.64 

115.2 

0.4305 

10.56 

0.14 

25.2 

0.3734 

9.160 

0.66 

118.8 

0.4187 

10,27 

0, 16 

28.8 

0.3953 

9.697 

0.68 

122.4 

0.4058 

9,954 

0.18 

32.4 

0,4145 

10. 17 

0.70 

126.0 

0.3919 

9.613 

0.20 

36.0 

0.4312 

10.58 

0.72 

129.6 

0.3768 

9.243 

0.22 

39,6 

0.4457 

10.93 

0.74 

133.2 

0.3605 

8.843 

0.24 

43.2 

0.4581 

11.24 

0,76 

136.8 

0.3429 

8.411 

0.26 

46.8 

0.4687 

11,50 

0.78 

140.4 

0.3239 

7.945 

0.28 

50.4 

0.4775 

11.71 

0.80 

144,0 

0.3036 

7.447 

0.30 

54.0 

0.4848 

11.89 

0.82 

147.6 

0.2817 

6.910 

0.32 

57.6 

0.4905 

12.13 

0.84 

151.2 

0.2582 

6.334 

0.34 

61.2 

0.4947 

12.12 

0.86 

154.8 

0,2330 

5.715 

0.36 

64.8 

0.4977 

12.21 

0.88 

158.4 

0.2060 

5.053 

0.38 

68.4 

0,4994 

12.25 

0.90 

162.0 

0.1771 

4.344 

0,40 

72.0 

0.5000 

12.27 

0.92 

165.6 

0.1461 

3.584 

0.42 

75.6 

0.4995 

12.25 

0.94 

169.2 

0.1131 

2.774 

0.44 

79.2 

0.4979 

12.21 

0.96 

172.8 

0.0778 

1.908 

0.46 

82.8 

0.4953 

12.15 

0,98 

176.4 

0.0401 

0.984 

0.48 

86.4 

0.4917 

12.06 

1.00 

180.0 

0.0000 

0.000 

0.50 

90.0 

0 ,4878 

11.97 

Model  4620 

Wetted  Surf.  Coeff.  = .7324 

Serial 

40050060 

-73 

LCB,  X = 

0.4456 

Formula : 

L/D  = 7.339 

y =aj^x+a2X  H-a^x 

+a,x  +a.x 
4 5 

+a.x 

D 

Model  Particulars: 

where : 

^1  = 

1,000000 

Length,  ft 

= 15.0000' 

^2  = 

1.137153 

Diameter,  ft 

Nose  radius,  ft 

= 2.044  (24.53  in) 

= 0.1392  (1.670  in) 

a-  = -10.774885 

Tail  radius,  ft 

= 0.0000 

3 

Wetted  Surf. , f t^ 

= 70.55 

a,  = +19.7»4286 

4 

Volume,  ft3 

= 29.53 

a = -16.792534 

LCB,  ft 

= 6.6840 

^6  = -^ 

5.645977 
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Landweber,  L.  and  M.  Gertler,  "Mathematical  Formulation  of  Bodies  of 
Revolution,"  DTMB  Report  719  (Sept  1950). 
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Figure  7a  - Stern  Appendage  Configuration  - NSRDC  Model  4620 


Figure  7 - Representation  of  NSRDC  Model  4620 
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Figure  8 - Lifting  Line  Circulation  and  T''^ch  Representing 
Loading  of  Propeller  3638 
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Figure  9 


Control  Points  and  Propeller  Induced  Velocities 
on  Appended  Series  58  Body 
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The  computed  values  of  thrust  deduction  and  the  measured  result  ob- 
tained from  model  resistance  and  self-propulsion  experiments  are  shown  in 
Table  3.  Since  no  correction  for  increased  frictional  drag  has  been  includ- 
ed, the  potential  flow  prediction  of  (1-t)  should  be  higher  than  the 
measured  value*.  The  agreement  between  the  two  calculated  values  (pressure 
integration  and  Lagally  theorem)  indicates  good  numerical  accuracy.  Also 
note  that  the  correction  to  body  singularity  strengths  changes  the  thrust 
deduction  fraction  by  less  than  3 percent.  The  comparison  between  lifting- 
line and  lifting-surface  predictions  shows  that  the  influence  of  blade 
thickness  and  chordwise  loading  reduces  the  thrust  deduction  fraction  by 
only  4 percent.  The  agreement  between  theory  and  experiment  is  quite  en- 
couraging and  is  considered  to  be  within  experimental  accuracy. 

It  is  believed  that  the  discrepancy  between  Beveridge's  result  and  the 
current  work  lies  in  the  derivation  of  the  propeller  sink-disk  strength. 
Beveridge  used  the  induced  axial  velocity  at  a lifting  line.  The  present 
method  (equation  27)  properly  averages  the  velocity  field  at  the  disk  plane, 
resulting  in  a lower  sink  strength,  as  shown  in  Figure  10. 

It  is  also  of  interest  to  examine  the  distribution  of  the  interaction 
force  over  the  afterbody.  This  is  a function  of  both  the  body  sectional- 
area  distribution  and  the  breadth  and  intensity  of  the  propeller  disturb- 
ance field.  The  longitudinal  distribution  of  the  thrust  deduction, 
(dt(x)/dx),  acting  on  the  body  only ^ is  shown  in  Figure  11.  The  signifi- 
cant contribution  is  over  the  last  25  percent  of  the  body  length  with  50 
percent  concentrated  in  the  last  6 percent  of  the  length. 

EXAMPLE  2:  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRESSURE 
DISTRIBUTIONS  ON  THREE  BODIES  OF  REVOLUTION 

Wind-tunnel  experiments  were  recently  conducted  to  determine  the  flow 

9 

characteristics  of  three  streamlined  bodies  of  revolution.  Measurements 
of  the  afterbody  pressure  distribution  with  and  without  a stern  propeller 
in  operation  were  obtained  for  comparison  with  analytical  predictions.  The 
profile  of  the  parent  body  (DTNSRDC  Model  5225-1)  is  shown  in  Figure  12 


*The  added  frictional  resistance  is  expected  to  increase  the  value  of  t 
by  no  more  than  5 percent,  based  on  Reference  9. 
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1.6 


Figure  10  - Radial  Distribution  of  Propeller  Sink  Strength 
Derived  From  Lifting  Line  Calculations 


Figure  11  - Calculated  Longitudinal  Distribution  of  Thrust 
Deduction  Force  on  Model  4620  (Body  Only) 


Figure  12  - Representation  of  DTNSRDC  Models  5225-1,  5225- 

and  5225-3 


Figure  12a  - Profiles  of  Axisymmetric  Afterbodies 


Figure  12  (Continued) 
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Figure  12  (Continued) 
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along  with  the  quadrilateral  representation  used  in  the  calculations.  The 
other  two  bodies,  DTNSRDC  Models  5225-2  and  5225-3,  differ  from  the  parent 
in  afterbody  shape  only  and  are  also  illustrated  in  Figure  12. 

One  7-bladed  propeller  model,  DTNSRDC  Model  4577,  was  used  for  the 

three  body  experiments.  The  principal  characteristics  and  geometry  of  the 

propeller  are  summarized  in  Table  4.  Pressure  measurements  were  obtained 

at  two  advance  ratios,  J = 1.25  and  J = 1.07.  Since  the  propeller  was  not 

operating  in  its  design  wake,  it  was  necessary  to  predict  the  load  distribu- 

21 

tion  using  an  Inverse  performance  calculation.  This  was  an  iterative 
process,  starting  with  the  measured  nominal  wake,  computing  a load  distribu- 
tion, recomputing  the  wake,  and  so  on.  Convergence  was  achieved  after  two 
iterations. 

The  computed  and  measured  differences  in  afterbody  pressure  distribu-  ' 

1 

tion  are  shown  in  Figures  13a,  13b,  and  13c.  The  agreement  is  excellent  1 

for  body  5225-1  and  very  good  for  body  5225-2.  There  is  a marked  dis-  | 

crepancy  for  5225-3.  This  body  was  intentionally  designed  to  have  boundary  ■ 

layer  separation  in  the  absence  of  the  propeller  with  the  hope  that  the 
propeller  influence  would  reattach  the  flow.  However,  flow  separation 
occurred  at  x/L  = 0.92  which  is  too  far  upstream.  The  thrust  deduction 
fraction,  was  derived  by  integrating  the  measured  pressure  difference 

and  dividing  by  the  calculated  propeller  thrust  (no  self-propulsion  experi- 
ments were  conducted).  The  results  are  compared  in  Figure  13  with  the 
theoretically  computed  values,  The  agreement  for  body  5225-3  is 

considered  fortuitous  in  view  of  the  difference  in  pressure  distribution. 

The  thrust  deduction  for  body  5225-2  is  larger  than  for  5225-1  which  is 
expected  on  physical  grounds  since  the  fuller  afterbody  places  the  frontal 
area  closer  to  the  propeller.  Body  5225-3,  which  has  the  fullest  after- 
body, is  actually  finer  in  the  immediate  vicinity  of  the  propeller;  hence, 
the  calculated  thrust  deduction  is  also  lower  than  for  body  5225-2. 


TABLE  4 - GEOMETRY  OF  PROPELLER  4577 


Number  of  Blades 

Expanded  Area  Ratio 

Section  Meanline 

Section  Thickness  Distribution 

Rake  Angle,  deg 

Skew,  deg 

7 

0.584 

NACA  a = 0.8 

BUSHIPS  Type  II^^ 

6.964 

30 

r/R 

C/D 

t/C 

P/D 

0.2106 

0.171 

0.235 

0.823 

0.0014 

0.3 

0.177 

0.209 

0.980 

0.0175 

0.4 

0.182 

0.182 

0. 151 

0.0288 

0.5 

0.185 

0.158 

0.243 

0.0337 

0.6 

0.185 

0.135 

1.264 

0.0341 

0.7 

0.180 

0.116 

1.248 

0.0311 

0.8 

0.164 

0.0995 

1.206 

0.0246 

0.9 

0.132 

0.0875 

1.157 

0.0141 

1.00 

0.069 

0.0813 

1.108 

0. 

0 / 

^Brockett,  T.,  "Minimum  Pressure  Envelopes  for  Modified  NACA-66  Sections 
with  NACA  a = .8  Camber  and  BUSHIPS  Type  I and  II  Sections,"  DTNSRDC  Report 
1780  (Feb  1966). 
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Figure  13  - Measured  and  Calculated  Pressure  Distributions  on 
DTNSRDC  Models  5225-1,  5225-2,  and  5225-3 


Figure  13b  - Body  5225-2 


5225-2 


i 
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EXAMPLE  3:  BODY  OF  REVOLUTION  WITH  STERN  APPENDAGES 


This  example  provides  a further  comparison  of  computed  and  measured 
thrust  deduction  and  exhibits  the  separate  contribution  of  stern  appendages. 
The  analysis  is  applied  to  an  appended  body  of  revolution,  represented  by 
DTNSRDC  Model  5224-1.  This  model,  with  appendages  removed,  is  a geosym  of 
the  wind-tunnel  model,  DTNSRDC  5225-1,  described  previously  in  Example  2. 

The  afterbody  profile  and  appendages  are  illustrated  in  Figure  14,  together 
with  the  quadrilateral  representation  used  in  the  calculations.  Note  that 
each  control  surface  is  represented  properly  (in  contrast  to  Example  1),  but 
the  forebody  is  again  approximated  by  the  afterbody  image. 

The  propeller,  DTNSRDC  Model  4567A,  (a  geosym  of  the  wind-tunnel  model 
4577),  was  designed  specifically  for  the  appended  body.  All  calculations 
are  therefore  based  on  the  design  loading  characteristics.  The  circulation 
and  hydrodynamic  pitch  distributions  are  shown  in  Figure  15. 

Calculated  propeller  induced  velocities  on  the  body  are  given  in  Figure 
16  showing  both  the  lifting-line  sink  disk  and  lifting-surface  contributions. 
The  results  derived  from  a simple  point  sink  located  at  the  shaft  centerline 
are  also  presented.  At  distances  beyond  one  propeller  diameter  the  veloci- 
ties derived  from  the  point  sink,  lifting-line,  and  lifting-surface  propeller 
representation  agree  to  within  10  percent.  Closer  to  the  propeller,  the 
lifting-line  velocities  are  about  25  percent  too  large,  while  the  point-sink 
velocities  are  as  much  as  50  percent  too  large. 

Calculated  and  measured  values  of  the  thrust  deduction  are  presented  in 
Table  5.  The  higher  lifting-line  velocities  result  in  a prediction  of  t 
which  is  21  percent  larger  than  that  derived  from  the  lifting  surface  calcu- 
lation. The  predicted  value  of  (1-t),  while  lower  than  the  measured  value, 
is  within  experimental  accuracy.  The  contribution  of  the  stern  appendages, 
as  shown  in  Table  5,  is  24  percent  of  the  total  thrust  deduction.  (Although 
not  discussed  in  Example  1,  a comparable  value  of  25  percent  is  found  in  that 
case).  Finally,  the  distribution  of  interaction  force  on  the  afterbody  is 
displayed  in  Figure  17,  shown  as  the  integrated  force  as  a function  of 
distance  from  the  stern.  The  significant  contribution  is  over  the  last  30 
percent  of  the  body  length  with  50  percent  of  the  force  concentrated  in  the 
last  5 percent  (i.e.,  within  a distance  of  one  propeller  diameter). 
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AXIAL  AND  RADIAL  INDUCED  VELOCITIES 


CONTROL  POINTS 
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J-\> 
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-0.02 

Figure  16  - Control  Points  and  Propeller  Induced  Velocities 
on  Appended  Body  of  Revolution  (DTNSRDC  5224-1) 
Propeller  4567A 
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TABLE  5 - COMPUTED*  AND  MEASURED  THRUST  DEDUCTION  FRACTION  FOR 
APPENDED  BODY  OF  REVOLUTION  (DTNSRDC  MODEL  5224-1,  PROPELLER  4567A) 


METHOD 

1-t 

Measured 

(Resistance 

and  Self  Propulsion  Experiment) 

0,925 

Computed* 

Pressure  Integration  (Equation  (11)) 

t„  , „ , = 0.070 

Body  Only 

Lifting 

Surface 

t,  ^ = 0.022 

Appendages 

0.908 

Lagally  Theorem  (Equation  (13)) 

0.904 

Lifting 

Line 

Pressure  Integration 

0.888 

*Without  Corrections  for  Increased  Frictional  Resistance 
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(Cumulative  From  Stern)  on  Body  of  Revolution 
(DTNSRDC  Model  5224-1,  Propeller  4567A) 
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EXAMPLE  4:  APPENDED  BODY  OF  REVOLUTION  WITH  FOUR  RAKED  PROPELLERS 


In  this  example  the  thrust  deduction  analysis  is  applied  to  compare 
four  different  propellers  fitted  to  the  same  hull.  The  hull  is  an  appended 
body  of  revolution,  represented  by  DTNSRDC  Model  5224-2.  This  model  is 
identical  to  DTNSRDC  5224-1  described  in  Example  3 except  for  a slight  modi- 
fication to  accommodate  a large  propeller  hub  (see  Figure  18). 

The  four  propellers  are  wake-adapted  designs  with  equal  diameters  and 

thrust  loading  coefficients.  A comparison  of  principal  characteristics  is 

given  in  Table  6 and  propeller  geometry  is  presented  in  Tables  7 through  10. 

The  radial  distributions  of  bound  circulation  and  hydrodynamic  pitch  shown 

18 

in  Figure  19  were  derived  from  lifting-line  calculations  using  Lerbs' 
criterion  for  optimum  loading.  The  propellers  are  distinguished  by  varying 
amounts  of  skew  with  corresponding  "skew  back"  along  the  respective  geo- 
metric pitch  helices.  As  a result,  the  blade  rake  differs  considerably 
among  the  four  propellers,  as  illustrated  in  Figure  20.  Physically,  it  was 
expected  that  these  differences  in  blade  axial  positions  would  strongly 
influence  the  thrust  deduction. 

Calculated  and  measured  values  of  thrust  deduction  are  compared  in 
Table  11.  Both  the  theoretical  and  experimental  results  exhibit  a pro- 
nounced decrease  in  the  interaction  force  with  increasing  propeller  rake. 
However,  the  computed  values  of  1-t  are  substantially  lower  than  measured 
in  each  case.  This  discrepancy  is  very  likely  due  to  difficulties  in  the 
experiment  owing  to  the  weight  of  the  propeller  models.*  Comparison  with 
the  result  of  Example  3 provides  evidence  that  the  measured  values  of  t are 
too  high. 

The  calculated  contributions  of  the  stern  appendages  to  the  thrust 
deduction  are  given  in  Table  11.  These  vary  from  19  to  24  percent  of  the 
total,  which  is  comparable  to  the  result  of  example  3.  In  Figure  20  the 
rake  distribution  of  propeller  4567A  is  shown  as  a dashed  line  and  may  be 
compared  with  propeller  4487.  Since  propeller  4487  is  further  aft  on  the 
body  and  has  a larger  diameter,  the  thrust  deduction  fraction  should  be 
lower.  This  is,  in  fact,  predicted  from  theory  (0.059  versus  0.09),  while 
the  experiments  yield  a higher  result  for  4487  (0.109  versus  0.09). 

*These  models  were  approximately  5 times  heavier  than  typical  models  due 
to  material  (bronze  instead  of  aluminum)  and  size. 
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Figure  18  - Stern  Alteration  on  Model  5224- 


TABLE  6 - PRINCIPAL  CHARACTERISTICS  OF  FOUR 
RAKED  PROPELLERS  ON  MODEL  5224-2 


Propeller 

4487 

4486 

4488 

4489 

Number  of  Blades 

7 

5 

3 

3 

Propeller  Diameter 

Hull  Diameter 

0.743 

0.743 

0.743 

0.743 

Expanded  Area  Ratio 

0.707 

0.707 

0.707 

0.709 

Blade  Thickness  Fraction 

0.036 

0.047 

0.077 

0.083 

P/  at  0.7R 

° X 

1.4975 

1.503 

1.519 

1.498 

Axial  Position, 

0.8637 

0.9652 

0.9678 

0.9703 

Skew,  Deg 

32 

43 

72 

120 

*Here  is  defined  as  the  longitudinal  distance  from  the  bow  to  a 
propeller  reference  plane  passing  through  the  midchord  of  the  blade 
root  section. 
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TABLE  7 - GEOMETRY  OF  PROPELLER  4486 

Number  of  Blades  5 
Expanded  Area  Ratio  0.707 
Section  Meanline  NACA  a = 0.8 

Section  Thickness  Distribution  NACA  66  (DTNSRDC  modified  nose  and 
tail) 

Skew,  Deg  43 


r/R 

c/D 

t/c 

P/D 

f /c 
m 

0.2 

0.219 

0.191 

1.319 

0.0725 

0.3 

0.250 

0.148 

1.366 

0.0569 

0.4 

0.280 

0.117 

1.408 

0.0413 

0.5 

0.310- 

0.095 

1 .445 

0.0296 

0.6 

0.320 

0.081 

1A77 

0.0227 

0.7 

0.320 

0.070 

1.503 

0.0171 

0.8 

0.300 

0.061 

1.517 

0.0134 

0.9 

0.250 

0.060 

1.509 

0.0110 

1.0 

0.000 

0.  - 

1.475 

0.  - 
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TABLE  8 - GEOMETRY  OF  PROPELLER  4487 

Number  of  Blades  7 
Expanded  Area  Ratio  0.707 
Section  Meanline  NACA  a = 0.8 

Section  Thickness  Distribution  NACA  66  (DTNSRDC  modified  nose  and 
tail) 

Skew,  Deg  43 


r/R 

c/D 

t/c 

P/D 

f /c 

m 

0.2 

0. 163 

0.210 

1.295 

0.0725 

0.3 

0.184 

0.152 

1.347 

0.0502 

0.4 

0.206 

0.120 

1.392 

0.0357 

0.5 

0.222 

0.095 

1.431 

0.0262 

0.6 

0.231 

0.080 

1.467 

0.0193 

0.7 

0.229 

0.070 

1.448 

0.0149 

0.8 

0.214 

0.063 

1.517 

0.0113 

0.9 

0.239 

0.059 

1.521 

0.0091 

1.0 

0.000 

- 

1.508 

- 

TABLE  9 - GEOMETRY  OF  PROPELLER  4488 

Number  of  Blades  3 
Expanded  Area  Ratio  0.707 
Section  Meanline  NACA  a = 0.8 

Section  Thickness  Distribution  NACA  66  (DTNSRDC  modified  nose  and 
tail) 

Skew,  Deg  72 


r/R 

c/D 

t/c 

P/D 

f /c 
m 

0.2 

0.353 

0.157 

mBM 

0.0792 

0.3 

0.425 

0.132 

mmm 

0.0578 

0.4 

0.481 

0.111 

■MH 

0.0438 

0.5 

0,520 

0.095 

0.0344 

0.6 

0.540 

0.080 

mmm 

0.0276 

0.7 

0.537 

0.070 

1.519 

0.0231 

0.8 

0.502 

0.063 

1.519 

0.0199 

0.9 

0.416 

0.060 

1.501 

0.0180 

1.0 

0.000 

1.458 

- 

TABLE  10  - GEOMETRY  OF  PROPELLER  4489 

Number  of  Blades  3 
Expanded  Area  Ratio  0.709 
Section  Meanline  NACA  a = 0.8 

Section  Thickness  Distribution  NACA  66  (DTNSRDC  modified  nose  and 
tail) 

Skew,  Deg  120 


r/R 

c/D 

t/c 

P/D 

f /c 
m 

0.2 

0.368 

0.200 

1.693 

0.0816 

0.3 

0.426 

0.152 

1.600 

0.0635 

0.4 

0.478 

0.119 

1.540 

0.0510 

0.5 

0.520 

0.098 

1.510 

0.0425 

0.6 

0.543 

0.082 

1.503 

0.0354 

0.7 

0.538 

0.069 

1.498 

0.0300 

0.8 

0.500 

0.061 

1.483 

0.0258 

0.9 

0.415 

0.058 

1.448 

0.0233 

1.0 

0.000 

- 

1.400 

- 

4486,  4487,  4488,  and  4489 


TABLE  11  - COMPUTED*  AND  MEASURED  VALUES  OF  THRUST  DEDUCTION 
FRACTION  FOR  PROPELLERS  4486,  4487,  4488,  AND  4489  ON 
DTNSRDC  MODEL  5224-2 


Propeller 

Experiment 

^^"*^^Theory 

^'Body 

^Appendages 

4486 

0.909 

0.949 

0.039 

0.012 

4487 

0.891 

0.941 

0.045 

0.014 

4488 

0.936 

0.966 

0.027 

0.007 

4489 

0.964 

0.974 

0.021 

0.005 

*Without  corrections  for  increased  frictional  resistance 
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By  considering  the  longitudinal  distribution  of  thrust  deduction  (body 
only)  shown  in  Figure  21,  it  may  be  concluded  that  virtually  all  of  the 
interaction  force  develops  on  the  last  20  percent  of  the  body.  Within  this 
region  the  effect  of  rake  is  manifested  as  a change  in  propeller  induced 
velocities.  The  total  axial  induced  velocity  on  the  afterbody  is  shown  in 
Figure  22.  The  separate  lifting  line  and  lifting  surface  contributions  are 
displayed  in  Figure  23,  showing  corrections  for  skew  and  rake  to  the  essen- 
tial difference. 

The  effect  of  blade  thickness  is  to  decelerate  the  flow  ahead  of  the 
propeller  and  hence,  reduce  the  interaction  force.  In  the  present  example, 
this  effect  is  more  pronounced  even  for  the  least  raked  propeller  (4487) 
because  of  the  relative  volume  of  each  blade  compared  to  the  loading.  For 
propellers  of  equal  blade  volume,  the  influence  of  thickness  would  diminish 
with  increasing  rake.  Such  is  not  the  case  here  because,  for  increased  rake, 
more  material  is  required  for  equal  strength.  Thus  propeller  4489,  while 
raked  the  most,  also  has  the  largest  thickness  contribution  because  of 
blade  volume  (for  example,  it  is  approximately  twice  the  volume  of  propeller 
4487). 

CONCLUSIONS  AND  RECOMMENDATIONS 

An  improved  theory  and  computer-aided  numerical  analysis  have  been 
developed  for  predicting  the  added  resistance  (thrust  deduction)  arising 
from  the  propeller-hull  interaction.  The  theory  is  formulated  in  terms  of 
the  diffracted  potential  flow  about  the  hull  in  the  presence  of  the 
propeller- induced  velocity  field.  The  propeller  representation  is  derived 
from  lifting-surface  theory  which  includes  the  effects  of  blade  thickness, 
skew,  rake,  and  chordwise  load  distribution.  It  is  shown  that  these  effects 
may  be  regarded  as  corrections  to  the  moderately- loaded  lifting-line  (sink 
disk)  approximation. 

Important  features  of  the  analysis  and  comparisons  with  experimental  i 

results  are  illustrated  by  application  to  several  stern  propeller-appended  j 

body-of-revolution  configurations.  Based  on  these  sample  calculations,  the 
following  conclusions  are  drawn: 
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1.  For  nonseparating  hull  forms,  potential  flow  theory  predicts  the 
change  in  afterbody  pressure  distribution  and  total  thrust  deduction  within 
experimental  accuracy. 

2.  Calculations  of  the  thrust  deduction  fraction  based  on  pressure 
integration  and  Lagally's  theorem  agree  to  within  4 percent. 

3.  The  influence  of  the  propeller  rapidly  decays  with  distance.  In 
all  cases,  the  thrust  deduction  develops  within  the  last  30  percent  of  the 
hull  length. 

4.  Corrections  to  hull  source  strengths  (interference  or  diffraction 
effect)  changes  the  thrust  deduction  fraction  by  less  than  3 percent. 

5.  Lifting-surface  corrections  reduce  the  calculated  thrust  deduction 
fraction  by  as  much  as  20  percent  for  conventional  propeller  geometries  (or 
equivalently,  a 2 to  3 percent  reduction  in  required  thrust) . This  correc- 
tion will  increase  for  more  highly  raked  propellers. 

6.  Conventional  cruciform  stern  appendages  contribute  up  to  25  percent 
of  the  thrust  deduction  fraction. 

7.  The  use  of  an  afterbody  image  to  represent  the  forebody  reduces  the 
computational  effort  without  loss  of  accuracy. 

In  view  of  these  findings,  the  following  extensions  and  applications 
are  recommended: 

1.  Further  comparisons  with  available  experimental  data  should  be 
conducted  to  aid  in  refining  the  analysis. 

2.  For  purposes  of  preliminary  design  calculations,  the  thrust  deduc- 
tion analysis  based  on  the  lifting-line  model,  as  described  in  the  report, 
should  be  incorporated  into  existing  propeller  design  computer  programs. 

In  order  to  properly  account  for  rake,  it  will  be  necessary  to  develop  a 
curved  lifting-line  representation  of  the  propeller. 

3.  A parametric  study  should  be  undertaken  to  systematically  examine 
the  role  of  hull  form  and  propeller  characteristics  (e.g.,  diameter,  axial 
location,  radial  distribution  of  loading,  and  rake). 

4.  The  method  should  be  extended  to  include  other  propeller  configura- 
tions — contrarotating,  tandem,  ducted,  and  twin-screw. 

5.  The  method  should  be  extended  to  apply  to  surface-ship  configura- 
tions . 
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APPENDIX  A 

THE  FORCE  GENERATED  ON  A BODY  IN  POTENTIAL  FLOW  BY  AN 
ISOLATED  SINGULARITY 

Consider  a body  immersed  in  a uniform  stream,  V , with  an  isolated 

00 

point  singularity  located  external  to  the  body  at  = r^.  The  flow  is 
assumed  to  be  steady,  irrotational,  and  incompressible.  If  the  body  is 
represented  by  a surface  distribution  of  source  singularities  a(r  ),  the 
the  force  exerted  on  the  body  is  given  by  (Equation  (13)  in  the  text) 


(i^„) 


Q'  B 


B' 


B 


where  V-(r„)  is  the  velocity  induced  by  the  singularity  at  the  body  surface. 

X B 

This  relation  may  be  derived  by  considering  the  control  volume  illustrated 
in  Figure  A-1. 


Figure  A-1  - Control  Volume  for  Analysis  of  Force  on  Body 
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The  pressure  forces  acting  on  the  control  volume  surface  are  equal  to  the 
time  rate  of  change  of  linear  momentum,  M,  of  the  fluid  within  the  control 
volume.  For  steady  incompressible  flow,  one  obtains 
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On  the  body  surface,  it  is  required  that  V*n=0.  The  force  acting  on 
the  body  is 


= / 


F = I pn  dS 


or,  from  Equation  A-1, 
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The  pressure  is  related  to  the  velocity  by  means  of  the  Bernoulli  theorem 
as 


p » - p_  V*V 
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and  so  the  force  may  be  expressed  entirely  in  terms  of  the  velocity  on  the 
boundaries  and  S , as 
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It  is  now  necessary  to  develop  expressions  for  the  velocity  at  large 
distances  (I'rl  »)  and  in  the  immediate  vicinity  of  the  singularity 
(r-tg  -*■  0).  For  convenience,  the  velocity  may  be  regarded  as  the  sum  of 
three  terms, 
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The  body  velocity,  V (^) , is  given  by 
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where 


VB(r)  ~ Vg(#g)  - sS-v^^  Vg(f^g)  + O (e^) 
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If  the  singularity  is  a point  source  of  strength  m,  then 


or 


F = - p 
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which  is  the  desired  result.  The  solution  for  the  case  of  a doublet  singu- 
larity may  also  be  written  in  this  form,  although  the  derivation  is  con- 
siderably more  involved. 
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APPENDIX  B 

SINK  DISK  REPRESENTATION  OF  A MODERATELY  LOADED  PROPELLER 

The  circumferential  mean  axial  and  radial  velocities  induced  by  a 
moderately  loaded  propeller  (outside  the  slipstream)  can  be  related  to  the 
flow  generated  by  a sink  disk.  Consider  the  velocity  potential  associated 
with  a distribution  of  sources,  0(r,<P),  located  at  the  propeller  plane 
X = 0,  0 < r < R, 
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If  the  source  strengths  are  independent  of  angular  position,  <p  , the  integral 
may  be  rewritten  as 
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Integrating  by  parts  yields 
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where  it  is  required  that  a(P)  = 0.  The  axial  and  radial  velocities  are 
then  given  by 
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These  results  are  equivalent  to  the  induced  velocities  of  the  propeller 
(Equation  (26)  in  text)  provided  that 
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It  follows  that 
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where  the  constant  of  integration  is  set  so  that  a(R)  = 0.  For  typical 
propeller  loading  distributions,  a(r)  is  negative  over  most  of  the  disk, 
corresponding  to  a sink  distribution. 
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